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We present a scaling study of lattice QCD with O(a) improved Wilson fermions and a chirally twisted mass 
term. In order to get precise results with a moderate computational effort, we have considered a system of physical 
size of 0.75 3 x 1.5 fm 4 with Schrodinger functional boundary conditions in the quenched approximation. Looking 
at meson observables in the pseudoscalar and vector channels, we find that O(a) improvement is effective and 
residual cutoff effects are fairly small. 



1. Introduction 

Twisted mass QCD (tmQCD) was introduced 
in [ |]J to solve the problem of unphysical fermion 
zero modes in lattice QCD with two degenerate 
flavours of Wilson fermions. The fermion lattice 
tmQCD action reads: 



a 4 V ij){x) (D + m + ip q <y 5 T 3 ) ip(x) , (1) 



where D is the 0(a) improved Wilson lattice reg- 
ularization^] of If) and tuq = l/2n — 4 and p q are 
two real quark mass parameters. The lattice sym- 
metries and power counting determine the coun- 
terterm structure for renormalization and 0(a) 
improvement of the infinite volume theory. The 
physical equivalence of tmQCD to QCD in the 
continuum limit has been argued in [ [IJ and will 
be further detailed in [ || . 

We study tmQCD in a four dimensional box 
with Schrodinger functional (SF) boundary con- 
ditions. We wish to check that after on-shell 0(a) 
improvement a la Symanzik [ || a few renormal- 
ized quantities related to meson physics approach 
the continuum limit with residual scaling viola- 
tions that are small and compatible with being 
0(a 2 ). Following ref. [ we keep the spatial size 

*Based on a poster by M. Delia Morte presented at the 
International Symposium on Lattice Field Theory, August 
2000, Bangalore, India. 
1 See ref. [ 0] for notation and details. 



of the box moderately small (L ~ 0.75 fm) and 
choose T = 2L. 

The transfer matrix of lattice tmQCD with ac- 
tion (|l|) and csw = can be constructed in close 
analogy to ref. [ [h]] and turns out to be self- 
adjoint and strictly positive for |k| < 1/6 [ [jj. 
Following refs. [ |[ [ll[], the Schrodinger functional 
in tmQCD can be conveniently defined as the in- 
tegral kernel of an integer power T / a of the trans- 
fer matrix [ |jj . It has an Euclidean representation 
given by 



Z[p',p',C';p,p,C} = 



a -S[C7>>] 



(2) 



fields 



and can hence be considered as a functional of 
the fields at the Euclidean times and T. The 
structure of the transfer matrix implies that the 
boundary conditions for gauge and quark fields 
are the same as in the standard framework (/i q = 
0). The action S[U,^,ip] in eq. (JsJ) is the sum 
of the usual plaquette pure gauge action and the 
quark action, which takes the same form as in 
eq. (Q), provided one adopts the notational con- 
ventions of Subsect. 4.2 of [ 0). 

The tmQCD Schrodinger functional is expected 
to be finite after the standard couplings and 
boundary renormalization [ [|. At order a 
new boundary counterterms proportional to a/x q 
arise [ ||] . Their effect beyond tree-level will be 
neglected in this study. 
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2. Correlation functions 

In order to motivate our choice of observ- 
ables, we recall some properties of renormalized 
tmQCD. For simplicity, we first consider the the- 
ory with boundary conditions that are periodic 
in space and as specified in [ ^] in time direc- 
tion. The renormalization scheme (R) of tmQCD 
can be chosen [ [| to be consistent with the Ward 
identities of the flavour chiral symmetries: 

d„(A K ); = 2m R (P R ) a + 8 a3 i m (Sn)° (3) 
d,(V R ); = -2 m s ab3 (P R ) b (4) 

where the above relations are to be understood 
as operator insertions^ in the correlation func- 
tions and e is fully antisymmetric with e 123 = 1. 
Note that our chirally twisted parameterization of 
QCD involves two renormalized mass parameters, 
m R and /iR . 

It can be shown [ |l|, || that, within a certain 
subset of the renormalization schemes that pre- 
serve the above Ward identities, the on-shell cor- 
relation functions of tmQCD can be mapped onto 
the ones of standard QCD in a related renormal- 
ization scheme. In the latter scheme the flavour 
chiral Ward identities (at the operator level) read: 

d^An^ = 2m' R (P R r , 9 M (V R )« = (5) 

for a — 1,2,3 with the mass mj satisfying: 

K) 2 = m K + fi K , tan(a) = ^ R /m R . (6) 

Here a is an unphysical angle which just specifies 
the quark mass parameterization. An example of 
this mapping is given by: 

((A R ) 2 (x)(P R ) 2 (y)) c {g R ,m' R ,0} = 
((A' R )l(x)(P R f(y)) c {g R ,m R , m }, (7) 

where i/i/ and by definition: 

= cos(a)(A R )«+ S m(a)e ab3 (V R )l, 
(PkT = (PrT a = 1,2. (8) 



It is understood that in eq. (0) the l.h.s. 
must be evaluated at renormalized couplings 



2 The quark bilinears are defined in the standard way, with 
flavour matrices r a /2 (for a = 1, 2, 3) or t° = 1. 



{g R , m' Rl 0} and the r.h.s at renormalized cou- 
plings {<7 R , m R , n R }, with the quark mass cou- 
plings related via eq. (||). Note that the above 
mapping takes the same form as at the classical 
level, where it is obtained via [ [jj 
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In the following we will also need to consider, 
again for a = 1, 2 only, the primed fields: 



{V R )l = cos(a)(V R )l + sin(a)e ab3 (A R ) b , 



(Tk) a 



a = 1,2 



(9) 



\1V V^^W/i 

Coming back to our scaling tests of tmQCD 
with SF boundary conditions, we focus on a few 
renormalized quantities that can be extracted 
from the following correlators: 



-({A' R )l(x)Ol)/Jfl 



(10) 



and (with sum over k 
1 



1, 2, 3 understood) 



1,22 



<{v R )l(x)Ql)l 



-^((T R ) 2 k0 ( X )Ql)/ 




(11) 



where and Q% are the SF boundary fields de- 
fined in ref. [ Q at zero Euclidean timeQ. In the 
above correlators the division by the square root 
of the quantity 

1 



Jl — 



L 6 



(o'tol) , 



(12) 



takes care of the boundary field renormalization. 

In the quantum mechanical representation of 
the renormalized correlation functions that are 
obtained from eq. (10) and eq. ([Tl"|), the only 
non-vanishing contributions (up to cutoff effects) 
arise when inserting states with vacuum quantum 
numbers at Euclidean times between xq and T, 
and states with pion and p-meson quantum num- 
bers, respectively, at Euclidean times between 
and xq. Let us consider for instance the SF cor- 
relator f R j±i ■ the matrix elements between inter- 



mediate states \k) and |n) 

(fc| cos{a)(A R )l - sm(a)(V R )l\n) 



(13) 



3 The analogous fields localized on the other SF time- 
boundary are labelled with a prime. 
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which enter its quantum mechanical representa- 
tion, are independent of the SF boundary condi- 
tions and also appear in the quantum mechanical 
representation of the r.h.s. of eq. (0). The equal- 
ity of the correlators in eq. (^|) implies the equality 
of their quantum mechanical representations for 
arbitrary values of (xo — yo)- From this equality 
it follows that in tmQCD the matrix elements of 
the form (|l3|) with the states |fc) and \n) carry- 
ing vacuum and pion quantum numbers, respec- 
tively, are non-zero. Analogous arguments hold 
for the remaining SF correlators in eqs. (|lO|)~(p"l|). 
The fields 0% and Qf. are chosen so that the cor- 
responding SF boundary states have a non-zero 
overlap with the pion and the p-meson states. 

3. Definition of the observables 

In this study we consider as observables some 
ratios of the previously introduced SF finite vol- 
ume correlators that are expected to approach a 
well defined continuum limin 



TOp S 


dofn 2 P ,(x ) 


Xq = 


--T/2, 


(14) 


tops 




X = 


--T/2, 


(15) 


my 


d Q k n 2 y,{x Q ) 
k R 2 v ,(x ) 


X Q = 


--T/2, 


(16) 


TOy 


d k R 2 T ,(x ) 


X = 


--T/2, 


(17) 


Vps 


= Cps/r 2 a'(iz:o) , 


x 


= T/2, 


(18) 


Vps 


= Cps/r,A'( 3; o) i 


xo 


= T/2, 


(19) 




= c v k R.y( x o) , 


x -- 


= T/2, 


(20) 




= Cvk R 2 Y > (x ) , 


x -- 


= T/2 . 


(21) 


As T 


= 2L — > oo, tops 


and 


my yield estima- 



tors of the pion and p-meson mass, respectively. 
The constants Cps and Cy are defined as in [ 
Q in terms of tops and TOy . Cps is such that 
77ps — » as T = 2L — > oo. The quantity r/y is 
not related to the decay of the p-meson because of 
its unphysical normalization. Analogous consid- 



erations hold for the alternative scaling quantities 
labelled with a "tilde". 

The /p- and fcp correlators entering our observ- 
ables are defined in eqs. (|lO|)— ( 12 ) via eqs. (||)-(||) 
and 



(A K )l = ZAil + bAam^A^, 

(V K )l = Z v (l + 6varn q )(Vi)« , 

(p R ,r = zp(i + fep fl TO q )(Pi) Q , 

(T R )J„ = ZTil + bTam^T^, 

where (restricting attention to a = 1, 2) 

(AOJ = A* + CAad^P a + af i q b A e ab3 V^ , 

{Vi)l = V^ + cyad^ + a^bye^Al 

(Pi) a = P a , 



(22) 



rpa 



c T a(8,V u a - d„V?) 



(23) 



For the definition of the bare lattice fields A?., 
P a , V" and T£ v we follow Sect. 2 of ref. [ |. 
The suffix I refers to the 0(a) improvement of 
these fields, which together with 0(a) improve- 
ment of the action, eq. (Q), implies that -in the 
limit T — > oo- the observables in eqs. (pT|)- (pl|) 
deviate from their continuum limit by 0(ar) cut- 
off effects. The b coefficients multiply countert- 
erms that are needed to subtract (bulk) cutoff 
effects of order ap q . 

4. Results 

We choose non-perturbative mass-independent 
renormalization conditions for the couplings and 
the observables of interest. As L/a is increased, 
the relation between j3 — 6/g 2 and a/r [ |l2| is 
employed in order to keep the physical size of the 
box fixed at L — 1.49ro, with ro being a hadron 
scale of order 0.5 fm. The mass parameter k is 
tuned so to keep fixed the PCAC renormalized 
mass in units of L: 



def L I Za < 
Ltor = - < — - 
a Zp 



(24) 



The symbol 9 M denotes the symmetric lattice derivative. 



where: 



= -{(A R )l(x)Ol) 
= -((PrY(x)OD 



(25) 
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The mass parameter n q is chosen so to fulfill the 
condition: 

L/x R d = Z^L^ = 0.153 . (26) 

This choice is justified by the exact lattice Ward 
identity d*V 2 = 2[i q P 1 , where V® is the one- 
point split vector current that is conserved at 
/iq = and <9* denotes the backward lattice 
derivative. The angle a, given by tan(a) = 
Mr/wr, is close to ir/2. 



p 






L/r 


Lrrvn 


6.0 


0.01 


0.134952 


1.490(6) 


0.0228(23) 


6.14 


0.00794 


0.135614 


1.486(7) 


0.0203(30) 


6.26 


0.00659 


0.135742 


1.495(7) 


0.0201(23) 


6.47 


0.00493 


0.135611 


1.488(7) 


0.0180(24) 



Table 1 

The bare and renormalized parameters in our sim- 
ulations. As for L/a and L/iR see the text. 



The renormalization factors Za, Zy and Zp 
were determined in [ |(| and [ [l3|, where Zp is 
given as a function of (3 at the scale Lq — 1.436ro. 
As for the improvement coefficients, we employ 
non-perturbative estimates of csw, ca [ || and 
&v [ Moreover, we use 1-loop estimates of cy, 
ct, &a, bp, &t [ f§ as well as of the SF-boundary 
coefficients c t and c t [ [llj . The coefficients 6a and 
by vanish at the tree level and are of order 10 -2 
at 1-loop level [ 9. We have hence varied their 
value in the analysis in the range —0.2 0.2 for 
&A and —0.1 -T- 0.1 for by, without observing any 
statistically significant variation in our results. 

Our simulation points are listed in table |l|: the 
four values of f3 correspond to L/a = 8, 10, 12, 16. 
The quoted error on Lmp is purely statistical, 
due to the negligible uncertainty on ZxjZp. As 
Zp = Zp(Lq) is known with a relative uncer- 
tainty of about 0.5%, the condition ( |26| ) is imple- 
mented with this precision by adjusting /i q . We 
correct for small mismatches with the condition 
eq. ( p4[ ) using numerical estimates of the Ltor- 
dependence of our observables obtained via some 
extra simulations at /3 — 6. We finally extrap- 
olate our data to the continuum limit assuming 
convergence with a rate ex a 2 . Our fits are shown 
in figures 0- 0, where only statistical errors are 




0.000 0.005 0.010 0.015 0.020 

(a/L) 2 

Figure 1. Scaling ofmpgL and fhpsL. 




2.20 

0.000 0.005 0.010 0.015 0.020 

(a/L) 2 

Figure 2. Scaling ofmyL and my L. 

displayed. The results are compatible with 0(a 2 ) 
deviations from the continuum limit. 

Because of the finite value of T ~ 1.5 fm, 
the observables ((l|)-(j2l|) still depend on the SF 
boundary action and fields. This induces residual 
O(a) effects due to the imperfect knowledge of Ct 
and Ct and a further coefficient^ that is associated 
with 0(a/j, q ) effects [ ||. Following [ we have 
performed a few extra simulations at /3 = 6 and 
f3 = 6.26 with values of c t — 1 and c t — 1 that are 
about 2 and 10 times, respectively, larger than 
the 1-loop values. The discrepancies with the 
previous results are at most about two standard 
deviations (at f3 — 6) and can be considered neg- 

5 Due to LfiR = 0.153 1, this coefficient was set to its 
tree level value and never varied. 
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ligible for the purposes of this study. 
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Figure 3. Scaling of rjpgL and r/pgL. 
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Figure 4. Scaling ofrjv and rjy. 

In table || we compare the estimated continuum 
limit value of our observables with the value at 
(3 = 6. The deviations from the continuum limit 
appear to be fairly small and not larger than the 
ones found in an analogous study at fi q = [ [7j] . 

5. Conclusions 

In the parameter region specified by (3 > 6, 
L/z R = 0.153 » Lm K = 0.020 and T = 
2L ~ 1.5 fm the O(a) improvement programme of 
tmQCD has been successfully implemented and 



mpsL 






Vv 


1.80(3) 


2.62(4) 


0.559(7)[1] 


0.164(5)[1] 


3.7% 


2.8% 


2.6% 


4.8% 


fhpsL 


m\[L 




Vv 


1.66(1) 


2.29(3) 


0.581(5)[1] 


0.200(6) [1] 


2.7% 


4.4% 


2.5% 


7.7% 


Table 2 



Continuum limits and deviations from (3 = 6. Er- 
rors due to Zx-uncertainties in square brackets. 

tested for a few meson observables. The renor- 
malization of the twisted mass parameter is easy 
in practice, as it can be traced back to the renor- 
malization of the non-singlet pseudoscalar density 
in the massless theory. 
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